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I. INTRODUCTION 


As more and more refinements have been developed on composite 
material during the past several years, the use of these materials in 
structural applications has become commonplace, if not the material of 
choice, in industry today. The reasons for this are well-known, the most 
common being the excellent strength-to-weight and stiffness-to-weight 
ratios associated with many composite materials. In addition, through 
micro-mechanical analysis, the engineer can design the material to optimize 
the physical properties (such as strength, thermal or hygro properties) for 
the functional requirements of the structure. As in any new developmental 
process, improvements to the product are always sought and the recent past 
has seen rapid growth in fiber improvements with respect to fiber strength 
characteristics. These new and improved fibers are then being used in the 
production of todays composite material. Associated with the rapid growth 
in fiber improvements is the general lack of an experience data base 
regarding the reliability and maintainability of the structure throughout its 
operational life. 

The Advanced Composites Laboratory at the Naval Postgraduate 
School, Monterey, CA, is addressing these questions with a systematic 
approach to composite reliability research. Unlike traditional design and 
manufacturing processes utilizing materials with well known characteristics 
(such as aluminum) where very large design "data bases" are known which 


often permits "design by experience", the incorporation of composite 


materials 1n design and manufacturing is still in its infancy, thus an 
equivalent "data base" does not exist. Design of a composite structure must 
be based on a sound reliability model of the salient physical failure process. 
Under tension, the failure process of a fiber composite 1s sequential. 
Internal failures start with breaking of weak fibers within the composite at 
relatively low applied loads. The original stresses carried by the broken 
fibers are transferred by the matrix binder to be shared by the neighboring 
unbroken fibers. This load sharing mechanism provides local redundancy 
thereby delaying global catastrophic failure. With the addition of external 
loads, more fiber failure sites are created leading to clustering of the 
failure sites which leads to ultimate catastrophic failure. This failure 
process is modeled in the Local Load Sharing model. Also, questions 
dealing with size effect (strength 1s inversely proportional to physical 
dimensions) of the composite structure must be addressed during the design 
process. This investigation 1s a part of the research effort to develop a 
probabilistic model of composite reliability which is capable of predicting 
actual large dimension structures from data produced in small dimension 
laboratory samples. 

Composite reliability is a function of both the strength and life of the 
constituent fiber and matrix material. To establish a probabilistic model of 
a composite's reliability, an adequate number of experiments must be 
conducted which deal with the strength and life of not only the composite 
but also of its parent constituent fibers. The objective of this study was to 
gather composite strength reliability test data to be used in exploring the 


viability of predicting composite reliability using a Local Load Sharing 


(LLS) model produced from previously obtained single fiber strength 
statistics. This was accomplished by testing composite strands to failure 


and comparing the data with the LLS reliability prediction. 


II. BACKGROUND 


A. FIBER REINFORCED COMPOSITE STRAND FAILURE 
MODEL 


Graphite composite materials usually consist of high-strength, high- 
stiffness fiber encased in a ductile matrix. The strength characteristics of 
the composite are dependent on the modulus of the fiber, the modulus of 
the matrix and the effectiveness of the bond between the matrix and the 
fibers at their interface. The failure model used in this study was first 
introduced by B. Walter Rosen [Ref. 1]. The model addresses the failure 
of a composite, consisting of a matrix stiffened by uniaxial oriented fibers 
when subjected to a uniaxial tensile load parallel to the fiber direction. It 
assumes that the fibers have a statistical distribution of flaws or 
imperfections which results in individual fiber failures in the composite at 
various stress levels. The composite itself will fail when the remaining 
unbroken fibers, at the cross-section with the most fiber failures (weakest 
cross-section), are unable to carry the applied load. Therefore, in this 
model, the composite failure is dominated by the tensile fracture of the 
fibers. In analysis of the statistical distribution of flaws in an individual 
fiber, it has been shown [Ref. 2 and 3] that the fiber segments can be 
modeled as a series of links with the fiber as a chain. The chain (fiber) 
will fail when the applied load is large enough to cause failure of the 


weakest link. The statistical strength distribution of the links can therefore 


be approximated by extreme value distributions which may be expressed as 
a two or more parameter Weibull distribution. Such distributions, when 
represented in a Weibull Probability of failure plot appear as linear or 
piecewise linear lines. 

If one were to combine many individual fibers into a bundle and a 
uniaxial load was applied to it, individual fibers would begin to fail as the 
load on each fiber exceeded that of its weakest link. This initial failure 
load is much less than that of the bundle failure load. Once an individual 
fiber has failed, the load that it was carrying is now distributed over the 
remaining fibers within the bundle, increasing the applied load on each 
individual fiber. The broken fiber no longer carries any of the load. If 
the combination of the applied and the assumed load on the individual 
fibers remaining is less than its weakest link strength, then the fiber will 
remain intact. If it is greater, then the fiber will fail and the same transfer 
of loading takes place. This process will continue as the applied load is 
increased and bundle failure will occur once every individual fiber has 
been broken. 

When the fiber bundle is encased in a matrix, the failure process of the 
newly formed strand is changed from that of the bundle. The model that 
describes the strand failure consists of parallel fibers in an otherwise 
homogeneous matrix. Its fiber statistical strength distribution is assumed to 
be the same as before. Figure 1, Appendix A shows the failure model 
presented by Rosen [Ref.1]. It should be noted that in this model, the 
extensional stresses in the matrix are neglected relative to those in the fiber 


and the shear strains in the fiber are neglected relative to those in the 


matrix. This is a reasonable assumption for fibers that are very strong and 
stiff relative to the matrix. As a uniaxial load is applied to the strand, the 
fibers with the weakest links begin to fail just as described in the bundle, 
but unlike the bundle, the matrix provides a unique load sharing 
characteristic that enables the internal broken fiber to continue carrying an 
applied load. In the vicinity of an internal failure, the uniaxial load that 
was carried by the fiber is transmitted by shear through the matrix to 
adjacent fibers. These adjacent fibers now carry a load increased by a load 
concentration factor (Kr). The internal failures result in shear stresses that 
locally may attain very high values. At the fiber end (point of failure), 
stress (c) goes to zero and the shear (t) goes to a maximum value. 
Therefore, a portion of the fiber near the fiber end is ineffective in 
resisting the applied load. This distance is known as the ineffective length 
(5). Rosen has shown [Ref. 1] that the ineffective length (normalized to the 


fiber diameter) of this model to be 
e 
1 2 ә 
]-v /: 2 S 
6 1 í = ont ф i 








5, 2 e Gm 
y 
f 


where O; is the fiber diameter 
Vf is the volume fraction of fiber in the composite 
Ef is the modulus of the fiber 
Gm is the shear modulus of the matrix 


p is the fraction of the undisturbed stress value below which 
the fibers shall be considered ineffective 


As the applied load increases, the number of statistically distributed fiber 
failures increase, producing an accumulation of ineffective fiber lengths. 
When a sufficient number of these ineffective fiber lengths combine in the 
vicinity of one cross-sectional area of the strand, it results in a weak 
surface. This enables the onset of matrix/fiber interface debonding and/or 
crack propagation through the matrix, ultimately resulting in the composite 
failure. In summary, the models treat fiber failures as the result of a 
statistically distributed flaw (weakest link). Composite strand failures are 
the result of a statistical accumulation of the fiber failures over a given 
cross-sectional area. Therefore, fiber strength is dependent on length; that 
is, longer chains (fibers) have a higher probability of having a weaker link 
(flaw) than a shorter chain. This agrees with experimental data [Ref. 4] 
that demonstrates that fiber strength is a monotonically decreasing function 


of fiber length. 


B. LOCAL LOAD SHARING 

The Local Load Sharing model used in this study was developed by 
D. G. Harlow and S. L. Phoenix [Ref. 2 and 3] and is shown in Figure 2, 
Appendix A. The fiber-matrix composite is viewed as a planar structure of 
n parallel fibers partitioned into a series of m sections, called bundles, each 
with n fibers. Conceptually, each bundle (defined as the metric) can be any 
desired length. Harlow and Phoenix designated the metric length to be the 
ineffective length 8. The length of the material.is then simply the metric 
length times m. In the model, the bundles are considered to be statistically 


independent and the strength of the composite material is that of its weakest 


bundle. In addition, it 15 assumed that the strength of the fiber elements are 


statistically distributed and can be modeled by the Weibull distribution 
FİL) - 1 - EXP | : D | for L20 


where a is the shape (slope) parameter 

p is the location parameter 

Lis the independent variable 
As shown in Figure 2, Appendix A, if the bundle load is x (per fiber) then 
the adjacent surviving fiber elements carry a new load Krx where 

Kr = 1 + r/2 p 

and r is defined as the number of consecutive failed fiber elements 
immediately adjacent to the surviving element (on both sides). Kr is called 
the load concentration factor. The probability distribution of the strength 
of a composite with Weibull distributed parent fibers can be computed 
numerically using a recursive relation for different fiber failure cluster 
(Kr) configuration. Based on this exact calculation, an approximate 


representation 1s available [Ref. 5] and is summarized as follows: 
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Kr = 14/2 where r = 1,2,3 
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with the variables H,y,(L} = Probability of failure of a composite of length 
m times the metric length 


WiL} = Probability of & - failure in the bundle 
FL) s Probability of failure of the mn fiber 
elements 
k = 


The number of adjacent fiber breaks in a 
bundle 


dis A constant that is a function of a 
and Ky, r, a, B and L are as previously defined. Finally, it is important to 
note that for a poorly made composite with inferior fiber/matrix adhesion 


and voids, the failure process will approach that of the global load sharing 


of a bundle rather than the local load sharing of a composite. 


III. METHOD OF TESTS 


A. DESCRIPTION OF TEST EQUIPMENT 
1. Mechanical System 

The strand tensile strength testing was conducted utilizing an 
Instron Universal Testing Instrument, model 4206, shown in Figure 3, 
Appendix A. It is comprised of two major systems: a crosshead drive and 
control system, which applies tensile or compressive loading to a specimen; 
and a highly sensitive load weighing system, which measures the load 
applied to a specimen. 

An Instron 2512 series, 1000 Kilogram load cell which was used to 
measure applied loads, was mounted in the moving crosshead which is 
operated by two vertical leadscrews within the loading frame. 

The test articles (composite strand samples physically resembling a 
pencil lead) were held by two Instron Modular Hydraulic Wedge Collet 
Grips (series 2742), one mounted on the load cell and the other on the base 
of the machine. One-half inch diameter collets were used in the wedge 
grips for mounting the specimen into the grip. The specimens were held 
stationary in the grips by hydraulic pressure obtained by a Instron Electric 
Pump Pressurization System (model A7154, Rev D, #5215) shown in 
Figure 3, Appendix A. The unit is self-contained with gripping pressure 
fully adjustable and independent of system pressure or fluctuations such 


that gripping force remains constant on the sample under test. 
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2. Control and Data Acquisition System 

The control and data acquisition of the testing instrument was 
automated with the use of Instron Series IX Materials Testing System. This 
system is a software package that interfaces the test instrument with a 
personal computer. The testing system components included the Series IX 
software, an IBM/PC-AT computer (1.5 MB of RAM, 30 MB hard disk) 
with an EGA monitor, an IEEE-488 (GPIB) interface, an IBM Proprinter 
II parallel printer and an HP 74704 plotter for graphic functions. 

The software, which is menu driven, presents the operator several 
options including computer controlled testing, reanalyzation of data, 
calibration of the test instrument, creating and or modifying the test 
method to be used in the computer controlled testing, and plotting of the 
raw data. Separate test methods were created for each gauge length of 


samples tested and are presented in Figures 4, 5, 6 and 7, Appendix A. 


B. TEST SPECIMENS 

1. Sample Composition 

The test specimens were produced from a Hercules Magnamite high 

Strength graphite, type AS-4. Both the fiber samples and the composite 
samples were fabricated from two specific production spools designated as 
008 and 019. Each composite strand consisted of 3000 fibers with a 
nominal cross-sectional area of .66 mm2. Samples with different gauge 
lengths were tested, the majority of which were 2 and 10 inches for spool 
008 samples, and 2 inches for 019 spool samples. Ten inch samples from 


the 019 spool were tested in previous investigations and are also included 
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for interpretation.. Two additional samples, one 5 inch and one 16 inch 
gauge length were also tested and this data used for compliance calculations 
in addition to normal data reduction. 
2. Test Method Development 

The samples were required to be made such that the strand could be 
gripped and held under load in the testing instrument. In addition, it was 
desired that the strand fail in the gauge length area and not in the gripping 
area. Two different testing methods were identified and investigated to 
accomplish these requirements. The first method involved securing the test 
strand ends in a copper tube (.0669 inch OD and .0390 inch ID) utilizing 
an epoxy type adhesive with the tubes separated by the desired gauge 
length. The two ends were then be placed in an adapter collet (from a 
manufactured jewler's lath having .0700 inch ID and .500 inch OD) which 
in turn were placed in the grips of the testing instrument. Several 
problems were encounter in the implementation of this method. First, the 
adapter collet/copper tube interface area was too small, producing a high 
stress concentration area on the strand itself upon gripping, causing failure 
of the strand in the grip. Also, the adapter collets were not designed to 
accept the type of loading the gripping procedure produced and they failed 
after a relatively low number of loadings. The second method involved 
using one-half inch aluminum bar stock to make 1/2 inch long pellets that 
were secured to the ends of the strands with the epoxy, resulting in a 
"dumbbell" type configuration of the test strand. Each pellet had a #60 
drill size hole (dimensions) machined into the center of it to accommodate 


the ends of the strands. This configuration allowed the test samples to be 


Iu» 


placed directly into the testing instrument, thereby simplifying the test 
procedures. Jigs, which were made out of aluminum angle iron, were 
designed and manufactured to maintain the pellets at the proper gauge 
length during the adhesive bounding process. At the same time, the jigs 
protected the samples from handling damage. This was the procedure that 
was followed to make the test articles and is explained in further detail in 
the following sections. 

Application of hydrostatic pressure, via the hydraulic grips, to the 
pellets produced a hydrostatic state of compressive stress in the 
sample/aluminim pellet interface. This increased the shear strength of that 
portion of the strand which faciliated the transfer of the external load to 
the fiber filaments within the free strand. Through trial and error, it was 
determined that a reduction of .0025 to .0035 inches of the outer diameter 
of the aluminum pellet upon application of gripping pressure was sufficient 
to produce an effective hydrostatic pressure. This amount of deformation 
was also the limit to prevent permanent deformation of the collets inside 
the grips themselves. To consistently obtain the desired yield 
characteristics, the pellets were all annealed at 450 deg Celsius for 1 1/2 
hours and then oven cooled overnight. This not only produced the desired 
deformation properties but also acted to rid the pellets of any residual oil 
and dirt deposited during the machining process. The annealed pellets 
were then tested by applying incremental steps of hydraulic pressure to the 
grips holding the pellet to determine the gauge pressure in the testing 
system that produced the desired deformation. This optimal pressure was 


estimated to be 500 psi indicated. 


3. Stress Wave Attenuation 

The intended purpose of a stress wave attenuator was to prevent 
secondary failures of the strand sample caused by the stress wave 
accompanied by the initial failure. Immediately after the initial failure of 
the sample the tensile energy stored in the unbroken portion of the sample 
Is released resulting in a tensile stress wave traveling toward the two ends 
of the grip. Upon reflection from the strand/pellet grip interface, the 
tensile wave is converted to a compressive wave. Since the strand sample is 
long and slender, the compressive wave almost always caused a secondary 
failure at the strand/pellet interface due to compress buckling. Additional 
secondary failure sites throughout the length of the sample were also 
fequently observed, therby precluded a definitive identification of the 
original initial failure site. It was experimentally observed that when the 
tensile stress was attenuated during the intial propagation stage, the 
secondary breakages could be eliminated. Because in the current testing 
configuration, the test sample was gripped in a vertical position, sand and 
oil were ruled out as candidates as stress wave attenuators, even though 
both of these have proven successful in the past.! The use of paraffin wax 
was also attempted but without success, as was modeling clay. The material 
adopted for these tests was a commercially available generic hot-melt 
adhesive that was applied, via a hot melt glue gun, to the entire gauge 


length of the test specimen. The general chemical constituents of the 


1 Wu/Nypiuk testing at Livermore Lab 
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adhesive include a styrene derivative, a hydrocarbon resin, artificial wax, 
and ethylene vinyl acetate coploymer resins. 
4. Test Specimen Preparation 

The composite strands were originally manufactured in 
approximately 18 inch sections and cut to lengths that were 1 1/4 inches 
longer that the desired gauge length by the use of a Dremel tool (using a 
cutoff grinding wheel attachment). The additional length was used to 
insure the strand extended through the entire length of the pellet. Once the 
required number of samples were cut (8 samples per batch), each end was 
cleaned by rinsing them through two separate acetone baths using a glass 
syringe. Additionally, the utensils (stirring stick and the syringe used for 
the placement of the epoxy) were also cleaned in the same manner. An 
epoxy that consisted of 55% by weight Dow DER-332 Epoxy Resin and 
45% by weight Texaco Jeffamine T403 was then made and stirred for a 
minimum of 3 minutes to ensure thorough mixing. The epoxy was then 
placed into a vacuum chamber to rid it of the air bubbles developed while 
mixing. This epoxy was used because of its particularly good wetting 
properties and its slow setting time. This allowed sample preparation to be 
conducted in a slow, methodical process reducing the possibility of 
damaging the samples during manufacturing. 

The pellets that were attached to the end of the strands (2 each) 
were then prepared by placing cellophane tape on three of the four ends of 
each set of pellets with small holes punctured through the tape in the center 
of each machined hole in the pellet. The use of a pin and microscope 


helped in placing the puncture hole in the desired location. This process 


PS 


was done to help center the strand in the pellet during the production of the 
test samples. The pellets were then placed in the jigs with the pellet with 
tape on both sides being the bottom pellet and the top pellet having the 
taped side down. They were held in proper place in the jig by rubber 
bands, allowing for the easy removal (cutting of the rubber bands) of the 
test sample from the jig at time of testing. The rear portion of the angle 
iron type jigs was machined away enabling access to the rear of the strand, 
making the application of the hot-melt adhesive easier. Once the jigs and 
the strands were prepared and ready, the epoxy was placed into a plastic, 
reusable syringe with a modified tip consisting of the bottom tip of a 5 mm 
mechanical pencil. This modification allowed the injection of epoxy 
directly into the top holes of the pellets until the epoxy flowed out the 
bottom. Since test samples were made in batches of eight, each of the 16 
pellets were filled with epoxy in this manner. The strands were then 
inserted into the pellets such that the strand extended through the entire 
length of each pellet. The jig was then placed in a vertical stand. Surface 
tension retaining the liquid epoxy in the cavity of the pellet alleviated the 
problem of inadvertently attaching the pellet to the jig. Once all eight 
samples were made, the stand was placed in an oven which was pre-heated 
to 58 deg Celsius (to accelerate cross-linking of the epoxy) and cured for 
24 hours. After curing of the adhesive in the pellet, each sample had the 
hot-melt adhesive applied, by hand, to the entire length of the strand. The 
samples were then ready for testing. Figure 8, Appendix A shows the 


completed samples. 
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C. TESTING PROCEDURES 
1. Alignment 

Before any testing was conducted, the test instrument's upper and 
lower grips were aligned by griping a precision ground 1/2 inch diameter 
steel drill rod in the grips with the grip surfaces separated by 
approximately 1/4 inch. A load of 1000 Kg (the load cell design limit) was 
then carefully applied and maintained (ensuring the load never exceeded 
design limits). To eliminate all of the mechanical play or looseness in the 
system, all of the rigid coupling attachments were then tightened under this 
tension,thereby setting the alignment of the grips. To ensure a consistent 
mechanical reference state for every sample tested, the grip alignment was 
checked and verified approximately one-half way through each testing 
interval. 

2. Calibration 

Load cell calibration was conducted before, during and after the 
testing period. Calibration was conducted using Instron calibration weights 
(traceable to National Bureau of Standards) with load increments of 5 
Kg. Calibration data is presented inTable I, Appendix B. In addition, a 
time drift or stability check of the system was conducted over a 24 hour 
period utilizing a constant load of 15 Kg. The greatest drift recorded 
during the period was within .17 per cent of the original load. As shown, 
the calibration remained well within the published specification limit of 
+1.0%. 
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3. Testing 

Testing of the prepared samples involved transferring the samples 
from the jigs to the test instrument and initializing and starting the test via 
the computer. The fully automated testing system recorded all the data and 
stopped the test upon failure of the sample. 

The placement of the sample into the test instrument required two 
people, one to assist the pellet of the strand in sliding into the collets 
internal to the grips and the other to operate the console that controls the 
crosshead displacement. The sample was first placed into the upper grip 
and inserted far enough into the collet such that the lip of the pellet was 
inserted beyond the lip of the collet. This ensured that the gripping force 
was distributed over the entire surface area of the pellet. A gripping force, 
measured by a hydraulic gauge pressure, of approximately 500 psi was 
then applied to the grip housing via the hydraulic pump. The crosshead 
was next lowered (manual mode, fast speed) such that the bottom pellet was 
approximately 1/4 inch above the upper surface of the lower grip housing. 
At this point, the speed select was changed to manual mode, slow speed and 
a load balance was conducted. Then the bottom pellet was inserted into the 
collet of the lower grip to the same relative position as previously 
described and gripping pressure applied. The crosshead was adjusted such 
that a slight negative load (0.0 to -.10 Kg) was indicated. The crosshead 
position reference was then reset to zero completing the installation. 

The computer was initialized, dependent on gauge length, utilizing 
the appropriate test method program. Sample identification numbers were 


developed such that they could be easily identified. For example, sample 
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01902-15 was the 15th sample of 2 inch gauge length of the AS4-019 class 
of composite strand that was tested. Once the required information of the 
interactive program was entered, the test was conducted. After testing to 
failure, the sample was removed from the test instrument and examined to 
determine the location of failure. A summary of all pertinent data 
collected during each test was then documented in the computer generated 
test report. A sample test report is presented in Figure 9, Appendix A. 
The raw data was finally stored on the hard disk and backed up on a floppy 
disk. 


IV. RESULTS 


A. TEST DATA 

1. Current Testing 

A total of 82 strands were tested to failure of which 56 were from the 
AS4-008 graphite spool and 26 from the AS4-019 spool. Load and 
displacement at failure data for all tests are presented in Tables II thru VI, 
Appendix B. Compliance calculation methodology for the correction 
factor applied to sample displacements due to system mechanical 
displacement is discussed in Appendix C. 

2. Previous Testing 

Data from previous tests of 008 and 019 composite strands that 
were part of the same physical set of the current test samples has been 
included in the data base for use in the analysis of data in this report. This 
was done to better define the reliability curve that the process of data 
reduction produces. The source of this data is from tests conducted at 
Lawrence Livermore National Laboratories by Dr. Edward M. Wu and 
Mr. Glenn Nypiuk [Ref. 6]. This data 1s presented in Tables VII and VIII, 
Appendix B. AS4-008 data from these tests are designated set À as 
compared to AS4-008 data set B which designates current test data. 
Fiber strength research for both (008 and 019) graphite spools was 

conducted by Dr. Edward M. Wu and Mr. Nypiuk [Ref. 6], Lt. David 
Keith Bell, USN [Ref. 7], LCDR Carl Engelbert, USN [Ref. 8] and Mr. Jim 


Nageotte of the Advanced Composites Laboratory, Naval Postgraduate 
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School. The fiber data from the above research, which is to be used in the 


final analysis of the present research, 1s contained in Ref. 9. 


B. HOT-MELT ADHESIVE CONTRIBUTION TO SAMPLE 
STRENGTH 


By using hot-melt adhesive as a stress wave attenuator, the question 
arises that, what contribution, if any, does the hot-melt adhesive make to 
the strength of the sample? This issue was resolved by testing the adhesive 
and determining its modulus which was then compared to that of the 
strand. Graphical representation of the hot-melt load-displacement data is 
presented in Figure 10, Appendix A. Knowing the cross-sectional area of 
the test sample, the modulus of the adhesive was calculated to be 2.845E3 
Ibs /in2. When compared to the modulus (approximated at 20E6 lbs / in2) 
of the composite strand, this 1s less than .014 per cent of the strand 
modulus. In addition to the modulus comparison, experiments were 
conducted to determine the effect the addition of the adhesive to the strand 
would make on the raw data. This was done by first testing a single strand 
(without adhesive) to a load of approximately 5040 of the failure load. 
This load application procedure was repeated several times to ensure that 
all the fibers that were weak and would break under a load up to and 
including the peak applied load, actually broke. This resulted in a fixed, 
repeatable, load-displacement curve for the test sample, indiscriminate of 
how many times a load (of the same peak magnitude) was applied. The 
hot-melt adhesive was then applied to the strand and the sample tested to 


the same load as before. As can be seen in Figure 11, Appendix A, the 
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load-deformation curve of the sample with the adhesive lies directly on top 
of the same curve of the sample without the adhesive. 

In analyzing the results of both of these experiments, it was 
concluded that in the context of the testing done to the composite strands on 
this project, that the hot-melt adhesive did not add to the strength of the 
composite strand and therefore the use of the hot-melt adhesive as a stress 


wave deflector did not effect the data. 


C. DATA ANALYSIS 
1. Data With Experimental Artifact 

In this type of data analysis, the data of those strands that failed 1n 
the gauge length portion of a test sample are considered the intrinsic load 
of that particular sample. The question then arises of how should the data 
of the samples that have multiple failure sites including one failure site at 
the grip/strand interface, be handled? We classified this data as censored 
data (due to the potential experimental artifacts which may have lowered 
the strength) . Based on this physical consideration, the failure load (E) 
for censored data samples must be alvvays less than and at best equal to its 
intrinsic value, that is, (L;) or Gay < (L;i). Our goal is to be able to use this 
censored data in combination with the intrinsic data in the data analysis. 
To achieve this, during data analysis a procedure using the Maximum 
Likelihood Estimator (MLE) was utilized. The MLE concept for a two 
parameter Weibull distribution is discussed in Appendix E. 

In data analysis, it is desired to plot the Weibull probability of 


failure of the material F* (Li) versus its intrinsic load. To do this, one 


must order the experimental data (failure load) in ascending order and 
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assign a rank to each point. In this analysis, expected ranking was used 


such that the rank of (L;) is 


až 
N+1 


fom 12.3... N 

where x is the numerical sequence of the failure load after ordering and N 
is the total number of samples tested. Because some of the data is censored 
and not intrinsic, a MLE analysis is conducted to estimate the model of, or 
the parameters of, the Weibull distribution associated with this particular 
set of data. A new probability of failure F(L;) is then determined for each 
point of censored data (Е) as shown in Figure 12, Appendix A. In 
analyzing F(L;), we note that if F (£) > F (L), it physically implies that the 
failure load (E) was greater than its intrinsic load (Lj) or (С) > (Li), which 
is inadmissible in the physical model. Therefore, in this case, for each 
point where F (Úi) =E (2), we assume (C) = (L;) and F(L;) = F(L;) or that 
the failure load is the intrinsic:load. If F(L;) <F (i) the value of F(L;) IS 
carried forward in the calculations. The new values of F(L;) are then 
combined with the original ranking of the intrinsic loads and the rank (with 
rank remaining keyed to the load) is reordered in ascending order. It 
should be noted that the MLE process should be performed on a sample set 
of a single dimension only. This is based on the fact that when the range of 
the probability of failure, F(L;) (ie, the sample size N} is small, the 
assumption of a linear Weibull distribution is a good approximation of the 
model. In other words, the model is assumed to be piecewise linear when 


the sample size is small. 


SE 


2. Data Considered Intrinsic 
This type of analysis assumes that the failure load of a particular 
test sample is the intrinsic load of the sample, or the sample failed at its 
statistically distributed intrinsic stress level or load. This assumption 
excludes any experimental artifact that may be associated with the data. 
This 1s the basis of the analysis used in this investigation. 
a. Data Interpretation 

Since both composite and fiber strength are functions of gauge 
length, a mathematical model must be determined to normalize all the data 
to one single metric size so that a direct comparison of all the data can be 
made. The standardization by a metric dimension can be chosen arbitrarily 
without affecting the relative relations. The metric selected by the current 
data analysis is one centimeter. The method for this standardization to a 
common metric dimension is discussed in Appendix D. The 
standardization used in the data analysis of this study is the vertical shift 
type. 

To determine if the data of like materials and gauge lengths 
are of the same set (identical test samples), the Weibull probability of 
failure plots are compared for each batch of tests completed. Figures 13, 
14, and 15, Appendix A show a comparison of each batch of data collected. 
In the analysis, batch one was the first half of the sample set which was 
compared to the second half of the sample set. Batch two compared all the 
odd number test samples to that of the even number test samples. The data 
band width is dependent on the experimental circumstances. A narrow 


band implies that the experimental techniques among the batches are 
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similar, justifying the merging of the different test batches into one data 
set. As shown, the data band width for each of the sample sets is narrow, 
indicating that the data belongs to the same set, which ¿s consistent to the 
actual experimental circumstances. By confirming that like batches are the 
same, they can be merged to produce single sample sets of like materials 
and gauge lengths. Figures 16 through 20, Appendix A presents the 
merged data of the sample sets as compared with the assumed Weibull 
distribution model obtained with the use of a Maximum Likelihood 
Estimator as described in Appendix E. The band width of the two sets of 
008 10" data (sets A and B) was then analyzed as shown in Figure 21, 
Appendix A and the two sets determined to be of the same set. The 008 
10" data was then merged and compared to the assumed Weibull 
distribution model as shown in Figure 22, Appendix A. The values of a 
and f produced by the MLE method for each set of data are presented in 
Table IX, Appendix B. All the data was then converted to a 1 cm metric 
format to enable a direct comparison with fiber data of the same metric. 
Figures 23 and 24, Appendix A present this data in terms of composite 
sample strength (per bundle) and Figures 25 and 26, Appendix A present it 
in terms of composite strength per fiber in a bundle. As shown in Figure 
27, Appendix A, in a direct comparison of composite strength (per fiber), 
the composite made with the 008 spool of AS4 fiber appears to be stronger 
than that made with the 019 spool of fiber even though both spools came 
from the same manufacturing line and were assumed to be identically 


produced. This is an indication that fiber spools that may have even been 
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produced in the same line, let alone the same lot, have to be regarded as 
having different statistical characteristics and treated as such. 

Much the same type of analysis can be done with the fiber 
data. Figures 28 and 29, Appendix A show the band width of the fiber 
data, which confirms that the data belong to the same set of like materials 
and gauge lengths. Figures 30 and 31, Appendix A present the merged 
fiber data as compared with the assumed Weibull distribution model. The 
values of the trimodal Weibull distribution for this data are presented in 
Table IX, Appendix B. As noted in Figure 30, Appendix A, the 019 fiber 
probability of failure is judged to be bimodal, or the lower tail having 
different Weibull distribution parameters than the upper tail. The modal 
shape of the data was determined by visual examination of merged 
individual ordered data plotted on a Weibull axis. The data was then 
identified, in this case, as either unimodal or bimodal. The fitting of the 
data to the upper and lower curves, including the identification of the 
lower transition point was performed heuristically. The MLE program 
was used to determine the values of e and B for both the lower and upper 
curves. The use of MLE is permitted for the bimodal case because it is 
considered to be a piecewise linear distribution in the Weibull plot with a 
distinct transition point. Figures 32 and 33, Appendix A present the 
Weibull probability of failure of the fiber when converted to a metric of 1 
cm. Figure 34, Appendix A shows the comparison of the fiber strengths 


and is consistent with Figure 27, Appendix A. ` 
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b. Local Load Sharing Analysis 

The data from Figures 25 and 32, Appendix A were combined 
to produce Figure 35, Appendix A. Likewise, Figures 26 and 33 were 
combined to produce Figure 36, Appendix A. These plots are used in 
conjunction with the Local Load Sharing (LLS) Prediction Model with 
extensions, to enable prediction of composite reliability from trimodal 
Weibull fiber statistics [Ref. 10]. The model requires two inputs, the fiber 
model parameters and the ineffective length, in order to predict the model 
for the composite reliability. The ineffective length, however, is not 
known, so the LLS program [Ref. 10] was written in a normalized (Xj/B) 
format. Therefore, when the model parameters of a particular set of fiber 
data are used in the program, the LLS curve that is produced is referenced 
to the normalized composite data of the same gauge length. In this 
analysis, one set of data (ie, the 008 data) was used in the model, since both 
the fiber and composite statistics (Figure 36, Appendix A) were known. 
Fiber statistics were entered using the model parameters in a trimodal 
Weibull distribution format. The model then generates a predicted value 
for the composite reliability from which the ineffective length is estimated. 
This is done by shifting the predicted reliability model horizontally to the 
right until the model approaches the value of the associated composite test 
data. Since the amount shifted is a function of the ineffective length of the 
composite, it is assumed to be the same for a composite made of the same 
fiber (AS4) but of different spools. The same magnitude of shift is then 
used in conjunction with the 019 fiber statistics and the LLS model to 


predict the 019 composite reliability. This prediction can be compared and 
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cross-plotted directly with Figure 35, Appendix A, making analysis of the 
LLS Prediction Model possible. The process can also be reversed by using 
the 008 data to obtain for the ineffective length and then predicting the 019 
composite reliability with the LLS model, which can then be compared 
with Figure 35, Appendix A. Figures 37 and 38, Appendix A presents the 
graphical comparison of the predicted composite reliability and the 
experimental test data for both 008 and 019 sample sets. 

As can be seen in Figures 37 and 38, the LLS composite 
reliability prediction exhibits outstanding correlation with the experimental 
data. In addition, it is sufficiently sensitive to predict the shifted composite 
strength due to the abnormally weak fiber tail associated with the 019 
fiber. It can also be inferred that the LLS reliability prediction provides 
the upper bound on reliability since the model assumes a well-made 
composite with a consistent ineffective length. This is further supported by 
the fact that the experimental data of the test specimens (which are 


considered well-made) correlated so closely with the predicted model. 


D. SUMMARY 

In summary, the major effort of this investigation was in the collection 
of a comprehensive set of composite data to be interpreted in conjunction 
with parent fiber data sets previously collected. The model for the data 
interpretation is the Local Load Sharing model which uses fiber statistics to 
predict the spatial density and clustering of the fiber failure sites. From 
the local stress concentrations that arise from load sharing, the probability 
of composite failure is predicted. This analytical prediction is then 


compared to the composite strength data sets obtained in this investigation. 
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Since the constituent fiber statistics are the benchmark for composite 
reliability prediction using the LLS prediction model, it is important to 
ensure that the experimental fiber data was obtained in such a manner that 
the data accurately reflects the statistical distribution of the fiber strength. 
This was accomplished in previous investigations by testing a sufficiently 
large sample population which ensures that the modal occurrence and the 
trimodal Weibull distribution parameters estimated are representative of 
the underlying fiber population. In addition, the use of multiple operators 
in conducting the tests also prevented the inadvertent skewing of the data. 
Finally different spools of fiber were used to identify normal 
manufacturing variables that are always present during the production 
process. 

To verify the prediction model with composite test data, the composite 
data must again be beyond question in proper representation of composite 
strength statistics. In this investigation, much time was spent in developing 
the test method procedures to ensure that the experimental data was 
accurate. The test sample gripping technique developed and utilized for the 
longitudinal tension testing was designed to prevent the normally dominant 
shear failure mode associated with tension testing from causing the test 
specimen failure. This was accomplished with a unique encapsulating grip 
that applyed hydrostatic pressure to the sample at the grip/sample interface 
which was induced by static compression. The use of the hot-melt adhesive 
to attenuate the stress wave that accompanies the sample failure, prevented 
the shattering of the sample and enabled the identification of the failure 


characteristics. In addition, different sample dimensions (gauge lengths) 


0 


were tested to extend the data range and enable the use of the methodology 
of data interpretation by shifting the weakest link. All this was successfully 
accomplished, producing great assurance in the experimental results. 

In interpretation of the data, for the case where the fiber is 
unimodal, the parameters for the two parameter Weibull distribution were 
estimated using the maximum likelihood estimator. For the case of data 
(019 spool of fiber) which displayed a bimodal distribution a model 
consisting of three piecewise linear Weibull distributions was used. The 
parameters for the trimodal Weibull distribution model were estimated 
from the ordered data set with censoring to the left and right respectively. 
The composite test samples that were used in this investigation were 
considered to be exceedingly well-made composite strands with uniform 
fiber and matrix adhesion and therefore, a uniform ineffective length. The 
ineffective length is controlled by the interfacial strength (a function of the 
manufacturing process) and the shear modulus of the matrix (materials). 
Composites with irregular fiber-matrix adhesion and voids have large 
ineffective lengths, the limiting case being a composite with no matrix (a 
bundle). It was observed that for this set of well-made composite, the 
Local Load Sharing Model of reliability prediction exhibited outstanding 
correlation with experimental data and was sufficiently sensitive to predict 
deficient composite strength due to a specific fiber spool with an 
abnormally weak lower tail. In addition, it provided an upper bound on 
the composite reliability. The documented weaker lower tail of the 019 
fiber was predicted in the composite reliability by the LLS model, 


verifying six orders of extrapolation of size in reliability prediction. It 
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should be noted that all analytical results were consistent with experimental 
observation, not only in trend but also in substantial agreement with 


magnitudes. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The results of this investigation suggest that the Local Load Sharing 
composite reliability prediction method using a trimodal Weibull 
distribution is a viable method of predicting composite structural reliability 
using fiber strength statistics. In design, this methodology provides a 
relation of strength data, based on small laboratory sample dimensions, to 
actual large structural dimensions for a fail safe design. In materials 
development, this methodology can be used to identify the parametric role 
of fiber and matrix strength properties for possible improvements. In 
manufacturing, this methodology can be used to prescreen fiber for zero- 
reject fabrications. In maintenance and repair, this methodology provides 
quantitative guidelines. 

The following areas were not fully treated in this investigation and are 
recommendations of follow-on research: 

1. That larger physical dimension test samples be tested such that 
next several orders of data extrapolation can be verified. 

2. That a mathematical formalism of the statistical methods to 
estimate the parameters for the trimodal Weibull distribution be 
established. 

3. "That a mathematical formalism of the statistical methods to 
estimate the best fit ineffective length based on the composite and fiber data 


be conducted. 
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4. That a closed form approximation of the Local Load Sharing 
composite reliability prediction method using trimodal Weibull distribution 


parameters be developed. 
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APPENDIX A 





Figure 1. Composite Failure Model 
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X - Load Per Fiber 

9 - Fiber Element Length 

Kr - Load Concentration Factor 
N - Number Of Fibers 

M - Number Of Bundles 


- Fiber Break 


Figure 2. Local Load Sharing Model 
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Figure 7. Test Method - 16" Gauge Length 
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Figure 10. Hot-Melt Adhesive Load Displacement Data 
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Figure 11. Load-Deformation Comparison of Hot-Melt 
Adhesive 
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Figure 16. 019 2" Weibull Model-Merged Composite Data Comparison 
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Figure 17. 019 10" Weibull Model-Composite Data Comparison 
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Figure 18. 008 2" Weibull Model-Merged Composite Data Comparison 
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Figure 19. 008 10" Set A Weibull Model-Merged Composite Data 


Comparison 
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008 10" Set B Weibull Model Data Comparison 
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Figure 20. 008 10" Set B VVeibull Model-Merged Composite Data 
Comparison 
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Figure 21. 008 10" Composite Band Width Comparison 
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008 All 10" Merged Weibull Model 


Data Comparison 
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Figure 22. 008 (All) 10" Weibull Model-Merged Composite Data 
Comparison 
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Figure 23. 019 Composite Data-Metric Relationship 
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Figure 24. 008 Merged Composite Data-Metric Relationship 
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Figure 25. 019 Composite Data (per fiber)-Metric Relationship 
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Figure 26. 008 Composite Data (per fiber)-Metric Relationship 
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008 - 019 Strength Relationship (Metric Size) 
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Figure 27. 008-019 Composite Strength Relationship 


60 
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Figure 28. 019 Fiber Data Band Width Comparison 
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Figure 29. 008 Fiber Data Band Width Comparison 
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019 Fiber Probability of Failure 
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Figure 30. 019 Weibull Model-Merged Fiber Data Comparison 
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Figure 31. 008 Weibull Model-Merged Fiber Data Comparison 
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Figure 32. 019 Fiber Data 


Metric Relationship 
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Figure 33. 008 Fiber Data-Metric Relationship 
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Figure 34. 008-019 Fiber Strength Relationship 
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019 Fiber/Composite Probability of Failure 
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Strength (per fiber) - [gm] 


Figure 35. 019 Fiber/Composite Probability of Failure 
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008 Fiber/Composite Probability of Failure 
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Figure 36. 008 Fiber/Composite Probability of Failure 
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Figure 37. 008 Composite Reliability Prediction 


Based on Fiber Statistics 
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Figure 38. 019 Composite Reliability Prediction 


Based on Fiber Statistics 
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APPENDIX B 


TABLE I. CALIBRATION DATA 


722 


Calibration Test Instrument Calibration 
Load (lbs) Load (Ibs) Drift (%) 
Before During After 
0 0.059 0.059 0.059 --- 
1102: 11.067 11.008 LT 127 +.54 
22.046 22.026 —17 228/35 +.27 
33.069 33.084 72 33.202 +.36 
44.092 44.092 43.974 44.092 +0.0 
55.116 55.101 55.041 55.160 +.11 
66.139 66.109 66.109 66.168 +.09 
71.102 7728807 7517 71:230 +.15 
88.185 88.185 88.125 88.244 +.07 
99.208 99.193 99.134 99.252 +.06 
110231 110.200 110.140 110.200 +0.0 


TABLE II. AS4-008 2 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


FAILURE LOAD EXPERIMENTAL! STRAIN AT 
(Kg) CERTAINTY FAILURE 
(mm) 
50.44 l .0330 
58.52 0 .0374 
56.54 l .0405 
51.76 0 .0410 
58.79 0 .0362 
55.09 0 .0355 
53.99 0 .0350 
53.88 0 .0351 
57.66 0 .0373 
55.84 0 .0357 
56.48 0 .0385 
58.71 l .0390 
57.61 0 .0362 
55.25 0 .0359 
55.30 l .0361 
58.25 0 .0388 
55.17 0 .0358 
57.53 0 .0375 
54.63 0 .0358 


Note 1: '0' - No known experimental artifacts 


]'- Known experimental artifact - intrinsic strength » failure strength 
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TABLE IL. 


FAILURE LOAD 
(Kg) 


56.00 
55.01 
53.66 
54.25 
55.28 


AS4-008 2 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


(CONT) 
EXPERIMENTAL! STRAIN AT 
CERTAINTY FAILURE 
0 0355 
0 .0352 
1 .0425 
0 0452 
0 .0367 


Note 1: '0' - No known experimental artifacts 


]' - Known experimental artifact - intrinsic strength > failure strength 
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TABLE III. 


FAILURE LOAD 


Note 1: 


(Kg) 


55.89 
23.99 
5577 
55.68 
57 
54.07 
54.98 
53.10 
5 
54.58 
47.49 
56.24 
55.60 
2157 
53.83 
54.17 
54.66 
35.57 
51.87 
52.81 
55.68 
54.44 
56.62 
54.17 
50 
35.96 
53.05 
52.56 
54.31 
53.15 


'0' - No known experimental artifacts 
]'- Known experimental artifact - intrinsic strength > failure strength 


AS4-008 10 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA (SET B) 


EXPERIMENTAL! 
CERTAINTY 


729 


— = Ço Ao = O = O A == OCC CO -»O» =>». 0». 00» =— 2. OOK CO 


STRAIN AT 


FAILURE 
(mm) 


.1710 
1799 
21:732 
.1856 
21759 
.1748 
217: 79 
.1702 
.1782 
.1770 
.1602 
.1775 
.1788 
.1658 
.1761 
.1776 
11 
.1181 
.1645 
.1710 
.1809 
11753 
.1794 
.1714 
1717 
21755 
21735 
.1688 
.1745 
.1763 


TABLE IV. . AS4-008 16 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


FAILURE LOAD EXPERIMENTAL STRAIN AT FAILURE 
(KG) CERTAINTY! (mm) 
53.42 0 2745 


Note 1: '0' - No known experimental artifacts 
] - Known experimental artifact - intrinsic strength > failure strength 
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TABLE V. AS4-008 5 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


FAILURE LOAD EXPERIMENTAL ! STRAIN AT FAILURE 
(Kg) CERTAINTY (mm) 
55.52 0 0896 


Note 1: '0' - No known experimental artifacts 


]'- Known experimental artifact - intrinsic strength > failure strength 


f 


TABLE VI. AS4-019 2 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


FAILURE LOAD EXPERIMENTAL STRAIN AT 
(KG) CERTAINTY! FAILURE 
(mm) 
54.50 0 0365 
57.74 0 0339 
58.79 0 0328 
60.24 0 0333 
58.93 0 .0372 
54.98 0 .0336 
51.60 0 .0341 
62.95 1 .0338 
47.25 0 .0320 
52.08 0 0344 
52.30 0 0362 
53.29 0 0361 
51.57 0 0361 
53.74 0 0345 
55.03 0 0355 
53.91 l .0347 
51.89 0 .0337 
54.58 0 .0403 
54.79 0 .0384 
51.03 0 .0332 
55.41 0 .0447 


Note 1: '0' - No known experimental artifacts 


] - Known experimental artifact - intrinsic strength > failure strength 
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TABLE VI. AS4-019 2 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


(CONT) 
FAILURE LOAD EXPERIMENTAL STRAIN AT 

(KG) CERTAINTY! FAILURE 
(mm) 

51.95 0 .0397 

54.76 1 .0448 

50.09 1 .0330 

51.46 0 .0349 

50.58 ] 0339 


Note 1: '0'- No known experimental artifacts 


]'- Known experimental artifact - intrinsic strength > failure strength 
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TABLE VII. AS4-019 10 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA 


INTRINSIC LOAD EXPERIMENTAL İl 
(Kg) CERTAINTY 


44.89 
42.38 
44.86 
47.15 
42.18 
50.60 
571 
49.08 
48.31 
43.47 
47.69 
45.61 
46.79 
50:35 
47.29 
48.44 
47.54 
45.18 
50.80 
44.48 
32:22 
46.53 
40.80 
48.05 


COD СУ (D OO 2 Q OO &2 QUO Cao]. oe) oS OUO ee 


Note 1: '0' - No known experimental artifacts 
] - Known experimental artifact - intrinsic strength > failure strength 
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TABLE VIII. AS4-008 10 INCH GAUGE LENGTH LOAD 
AND DISPLACEMENT FAILURE DATA (SET A) 


INTRINSIC LOAD EXPERIMENTAL | 
(Kg) CERTAINTY 


55.19 
45.85 
57:52 
43.24 
49.40 
53.66 
54.03 
54.28 
43.71 
1 
5170 
50.67 
48.65 
48.62 
54.70 
50.00 
51.34 
51:05 
59152 
51.47 
45.82 
44.04 
45.63 
48.25 
48.74 


cc a cC c O O C (О c C occ o orc clc Coo o o 


Note 1: "0" - No known experimental artifacts 
]'- Known experimental artifact - intrinsic strength > failure strength 


Š ] 


TABLE IX. WEIBULL DISTRIBUTION PARAMETERS 


Test Sample o D AL TLM AM BM TMU AU! 
008 Fiber(1cm) -- -- 5.23 3.75 -5.232476 3200 SS 
019 Fiber(1cm) -- -- 3.3] 9.62 5:08 23.75 3000 5.08 
008 Strand 2207777: = = 2s = si = 
(27) 

008 Strand II STE SM = së = T Gs H 
10" (Set A) 

008 Strand 26.00 56.09 -- == - Ss = = 
10" (Set B) 

008 Strand 155054 5 = = m x = = 
(Merged) 

019 Strand MS 502 SS es A = = = 
(2") 

019 Strand 17.48 48.34 = = £ = == ə 
(107) 

Note 1: Strand distribution is presented in a unimodal format. 


Fiber distribution is presented in a trimodal format to 
facilitate data entry into the LLS program. Abbreviations are 
A-Q B-B 

M - Middle 


T - transition point 


L - Lower U - Upper 
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APPENDIX C 


I. COMPLIANCE CALCULATIONS 


To obtain the most accurate load versus displacement representation 
possible, the raw cross head displacement data must be adjusted to take into 
account the compliance of the mechanical system itself. The displacement 
recorded (9) during data collection is a combination of the compliance of 


the mechanical system (0') and that of the test strand (9") such that 
d=d0'+0"' (1) 
Knowing 
du 2) 
where P is the load and k' is the equivalent spring constant of the 
mechanical system (which is a constant). For the strand itself, 


P o" .. 
os Ee or es or d''=— (3) 
pw £ 


where A is the cross-sectional area of the strand 
E 1s Youngs Modulus 
l is the gauge length of the sample 


Note that since the test strands are considered identical, 1/AE is a 


constant. 


8 3 


By substituting equations (2) and (3) into (1) for two samples of different 
gauge lengths, the resulting system spring constant is 
lz- l 


iie 4 
s dl, д 21, 


11. 


Since in data interpretation it's desirable to express k' as a function of a 
load (P), we take the case where the load is the same for each data point, 


thus P1 = P2 = P, therefore equation (4) simplifies to 


LAB) 
| d ə 
Knovving 
AP 
O) 


we can solve for d' at any point by integrating, thus resulting in 
uf 
A l 
ane f OY 
0 


By using equation (5), the system spring constant (k ') and compliance 
(6), both functions of load, were calculated at each point that a common 
value of the load occurred during the testing of two separate gauge length 
samples. Ideally, compliance data should be taken with a constant load 


control system instead of a constant cross-head rate control system for 
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direct data reduction. This feedback control configuration is not available 
for the testing system used in this investigation; as an alternative, the data 
from two different gauge length samples were manually reduced and 
common loads identified in order to conduct the analysis. Another 
approach that could be taken would be to curve fit the test data, establish an 
equation (displacement as a function of load) for each set of data, and then 
simply calculate the displacements at the desired loads. 

Four samples of varying gauge lengths (2, 5, 10, and 16 inches) were 
tested and analyzed to obtain the system compliance. The samples were 
produced as described in the Background section using AS4-008 composite 
strands. All combinations of 0'(P) computed from the four samples were 
compared and the resulting graphical comparison is presented in Figure 1. 
As shown, system compliance closely correlates in three (2-5 inch, 2-10 
inch, and 2-16 inch) of the six combinations. Since the calculated value of 
9'(P) at low loads for the 2-5 inch gauge length compliance is negative, the 
final compliance calculation excludes this result. Thus, the final system 
compliance was derived by averaging the 2 and 16 inch gauge length 
compliance data with the 2 and 10 inch gauge length compliance data which 
Is presented in Figure 2. The equation for the system compliance was 
obtained by using a third order polynomial curve fit through the averaged 


data and 1s 


Q'(P) 2 9.3062E-5 + 2.8145E-5 P + 7.9643E-8 P2 4- 4.7468E-10 P? 


It should be noted that this equation is valid only for loads up to 145 


lbs since extrapolation of the curve fit polynomials are not 


Š 5 


necessarily valid beyond this point. However, this is sufficient since 


the maximum loads observed during testing do not exceed this limit. 


0.04 


2-16 Data 
2-10 Data 
10-16 Data 
0.03 02-05 Data 
* 05-16 Data 
05-10 Data 


0.02 


0.01 


MECHANICAL SYSTEM COMPLIANCE (9!) 
[in] 





-0.01 
0 50 100 150 
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Figure 1. System Compliance Comparison 
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AVG MECHANICAL SYSTEM COMPLIANCE (9) 


{in] 


0.008 
0.006 ° AVG COMPLIANCE 
0.004 
0.002 
0.000 €” 


-0.002 
0 50 100 


LOAD [LBS] 


Figure 2. Avg. Mechanical System Compliance 
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APPENDIX D 


I. MERGING OF STRENGTH DATA FROM DIFFERENT 
SAMPLE DIMENSIONS 


There are two different mathematical normalizations of the strength 
data measured from samples of different dimensions (in the context of this 
investigation, different gauge lengths). The underlying assumption of the 
mathematical normalization is that the failure is a extreme value process 
(weakest link) and that the strength of each segment of the sample (the link) 
is Independent and Identically distributed. 

When data are represented in a transformed Weibull cumulative 
probability of failure, a horizontal shift of the data can be made which, in 
effect, shifts the location parameter of the assumed two parameter Weibull 
distribution model. This results in a shift of the lower tail of the Weibull 
probability curve down. The derivation for this correction of gauge length 
is presented in Appendix C of Ref. 7. It should be noted that this method 
can only be utilized on a two parameter Weibull distributed function 
(which 1s linear in the Weibull probability of failure coordinates) or a 
three parameter function that can be linearized into a two parameter 
function. 

A vertical shift of the data can also be made which maintains the same 
location parameter but shifts the Weibull probability curve up or down on 


its vertical scale in order to account for the change in probability associated 


8 8 


with a change in the metric dimension. It can be shown that the reliability 
of a chain (fiber) is equal to the reliability of a link in the chain (the metric 
dimension) raised to the power of the number of links in the chain (n), or 


in terms of probability of failures 


(1 - Fn(L)) = (1 - F(L))? 
where  Fpn(L) is the probability of failure of the chain and 
F(L) is the probability of failure of a link in the chain 


This assumes that the links of the chain are independent and identically 


distributed. The equation can be written as 


In(1 - Fn(L)) =n In (1 - F(L)) or 
In(-In(1 - Fn(L))) = In n + In(-In (1 - F(L))) 


which ıs observed to be in og form of x' =c + x. Therefore, the 
distribution of the chain preserves any and all distributions of the link and 
the distribution of the chain is simply a vertical shift of the link by 1n n (the 
physical size of the chain). Note that the use of In(-In(1 - Fn(L))) implies a 
series of weakest links (Type III, extreme value statistics) and if it 1s linear, 
then it implies a Weibull distribution. The same procedure used in the 
above mathematical derivation 1s also valid for a composite with load 


sharing. Applying these procedures results in 


In(-In(1 - Hmn(L2)) 7 In(imn) * In(-In(1 - W(L))) 
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where Hmn(L) is the probability of failure of a composite of length / 
(given mn) where 7 is the gauge length 
W(L) is the probability of failure of a composite at the 
desired metric length 
n is the number of fibers in the bundle 
m is the number of bundles (segments of the length 
of the metric) in the gauge length 


Note that this again is in the form of x" € c +x and that the transformation 


simply involves a shift in the vertical axis. Rearranging, we have 


In(-In(1 - W(L))) 2 InC-In(1 - Hmn(L))) - In@mn) 


This relation is used to transform strength data measured from different 
gauge lengths into one common gauge length. This normalization of all 
data into one reference metric dimension allows the merging or pooling of 
the aforementioned strength data for comparison with the analytical Local 


Load Sharing probability predictions. 
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APPENDIX E 


I MAXIMUM LIKELIHOOD ESTIMATOR 


A. BACKGROUND 

The Maximum Likelihood Estimator (MLE) method 15 а method of 
estimating the parameters of a model (curve) for a given set of data. There 
are several advantages of the MLE method over the least squares method. 
The least squares method assumes equal weighting of the data in the 
linearized domain and does not account for data clustering. In addition, it 
can only be used with exact (x = x; ) data, where x is the true realization of 
the data and x; is the experimental data. The MLE method, on the other 
hand, correctly weights the data by the probability of occurrence and it is 
applicable for exact (x = xj ), interval (xj < x < xj4 1), or censored 
(x > x; ) data. A data set may consist of one or more of the three classes 
of data. The MLE method is statistically sufficient, asymptotically 
consistent, asymptotically efficient , and converges to the underlying 
parameters. 

The principle behind the MLE method is to determine the model 
parameters of the assumed model that has the maximum likelihood, or 
probability, of repeating a particular set of data from all possible 
repetitions. In developing the mathematical representation of the MLE 


method, it is assumed that the combined probability of all the data 
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occurring during testing is the product of the probability of each data point 


in the set, or 


n 


TI f(x:,0) - ES ə) Ka) D 9) | M ,0)| (1) 


i= 1 
where f(x; ‚8) is the probability of one of n data points given the model 


parameters 0. This assumes that the data were collected from identical, 
independent samples. The likelihood (L) is therefore defined as this 


combined probability, or 


n 
L= II f(x:,0) (2) 


1-1 


One notes that to solve for the likelihood (L), the probability function 


f(x;,8) must be known, which includes both the data and the model 


parameters. The problem, of course, is that the parameters are not known. 

One approach to solving this problem is to conduct an iterative 
calculation process such that you assume values of o and p, and solve for 
(L). The likelihood (L) is then plotted as a function of « and p, resulting 
in a 3-D plot. Different values for e and B are then chosen and another L 
Is calculated and plotted. This process is continued until the plot can be 
analyzed and the peak or maximum likelihood determine. These values of 
o and B become the estimated parameters of the model. 


There are several problems associated with this process. First, since 


f(x;,0) « 1 , it implies that 
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n 
II f(x; ,8) | 
1-1 
and for large sample sizes, the computer's numerical capacity can be 
underflowed. In addition, the evaluation of where the maximum likelihood 
occurs involves the use of a concept dealing with Confidence Regions, 
which turn out to be coupled with associated Confidence Intervals for the 
parameters. The evaluation of (L) can then become quite complicated. 
Another solution to the problem is to solve directly for « and B. This 

solution maximizes the likelihood (L) analytically by taking the derivative 
(with respect to the model parameters) of equation (2) and setting them 
equal to zero. For the two parameter Weibull distribution, these 
mathematical steps can be performed explicitly resulting in two equations 
with the two unknown model parameters (x and B) These two equations 
are then solved simultaneously for the model parameters. As it turns out, 
the equations for œ and B are decoupled such that one parameter (normally 


a) can be optimized and then the other parameter (f) calculated. 


B. CENSORED DATA WITH KNOWN BUT DIFFERENT 
LOWER BOUNDS (WEIBULL DISTRIBUTION) 
1. Definition of Variables 


n z Total number of tests 


m = Number of tests of known realized strength (no uncertainty) 
x; = Exact data 

x; = Realized data where X; < x; (uncertain data) 

09 = Parameters of the assumed model 


OÜ) 


The likelihood is the product of all the individual probabilities, but the 
exact value for the uncertain data is unknown. What is known is that the 
exact data is at least as great as or equal to the uncertain data. Therefore 


the likelihood 1s 


L- А fix,0) = TE 2 f(«.8)dx]. ак 


CSL 


Or 


Knowing that 


Í dE yl J teseyax=1- FG) 


2 


the likelihood is then 


ER n 
L-|| n $x,e)|| n (1-F@,0)) 


i=l JEN 


For a 2 parameter Weibull distribution we know 


a-1 a 
OLX X 
f(x,9) = во) is | 
a D 





and 
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ğ xj 
] - F(x,,9) = exp à - = | 
B 


Therefore 





Or 
7: est? (qey 
L=|apB I dəl əfi II exp É 
E p : B 


Maximizing L with Respect to a and B 


2. 
To remove the multiplicative products from the above equation, 


we take the natural log (In) of L such that In(L) 2 L. Therefore 


L =m In(o) - ma In(B) + S (a - D) In(xj - вә x" B Y x; 


i=l 1-1 i=m 


To maximize L (equivalent to maximizing L) with respect to a and B , the 


partial derivatives are taken and set equal to zero, or 


ə and 0 


ÒQ op 


This results in 


oL Boe Y — “ə > ül 0] Elmo] | 
1-3 - 


Oo i=1 
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VIX B zə) 9 


and 





p= ub ə Y d (4) 


i=] i=l ] = 
Q= m In(x;) 
] m n ACC m i=l 
Y AX 
m ¡=1 i-m 


and 


96 


We note that the shape and scale parameters a and f are 
decoupled but a is implicit. Either an iteration calculation or numerical 
methods are used to solve for the shape parameter a. Once a is found, p is 
calculated directly. We further note that the equations for a and f can 


also be used for non-censored data by simply letting m = n. 
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